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Abstract 

The nonlinear evolution of two fluid interfacial structures like bubbles and spikes arising 
due to the combined action of Rayleigh- Taylor and Kelvin-Helmholtz instability or due to that 
of Richtmyer-Meshkov and Kelvin-Helmholtz instability resulting from oblique shock is inves- 
tigated. Using Layzer's model analytic expressions for the asymptotic value of the combined 
growth rate are obtained in both cases for spikes and bubbles. However, if the overlying fluid 
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is of lower density the interface perturbation behaves in different ways. Depending on the 
magnitude of the velocity shear associated with Kelvin-Helmholtz instability both the bubble 
and spike amplitude may simultaneously grow monotonically (instability) or oscillate with time 
or it may so happen that while this spike steepens the bubble tends to undulate. In case of 
an oblique shock which causes combined action of Richtmyer-Meshkov instability arising due 
to the normal component of the shock and Kelvin Helmholtz instability through creation of 
velocity shear at the two fluid interface due to its parallel component, the instability growth 
rate-instead of behaving as 1/t as t — > oo for normal shock, tends asymptotically to a spike 
peak height growth velocity ~ j^j^j^ j (An) 2 where At; is the velocity shear and At is the 
Atwood number. Implication of such result in connection with generation of spiky fluid jets in 
astrophysical context is discussed. 
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I. INTRODUCTION 

Ray leigh- Taylor (RTI) and Kelvin-Helmholtz (KHI) instabilities are associated with the perturba- 
tion of the interface of two fluids of different densities subject to the action of continuously acting 
acceleration (with respect to time) and under the action of velocity shear,respectively. The perturba- 
tion and the consequent instability may also be induced by a shock generated impulsive acceleration 
known as Richtmyer-Meshkov (RMI) instability. Such interfacial hydrodynamic instabilities occur 
in a wide range of physical phenomenon from those associated with problems on wave generation 
by wind blowing over water surface to problems related to Inertial Confinement Fusion (ICF) or 
astrophysical problems like that of supernova explosion remnant which belong to the domain of high 
energy density (HED) physics ^ . In ICF experiment HED plasmas may be created due to multi 
kilo Joule laser with a pressure ~ Mbar. In ICF, in addition to RT and RM instabilities nonspherical 
implosion generate shear flows; the later is also formed when shocks pass through irregular fluid 
interfaces. The KHI and shear flow effects in general are also of practical importance in a number of 
HED system. They should be considered in a multi shock implosion schemes for direct drive capsule 
for ICF, since KHI may accelerate the growth of turbulent mixing layer at the interface between 
the ablator and solid deuterium-tritium nuclear fuel. In HED and astrophysical system, it has been 
seen that structures driven by shear flow appear on the high density spikes produced by R-T and 
R-M instabilities^. They may develop in course of evolution of these instabilities El— El anc l cover 
enormous range of spatial scales from 10 17 cm for jets from young stellar objects to 10 24 cm for jets 
from quasars or active galactic nuclei ^ . Examples are suggested to be provided by pillars ("elephant 
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trunk" ) of Eagle Nebula which are identified with spikes of a heavy fluid penetrating a light fluid®, 
Another example in astrophysics is the Herbig-Haro (HH) object like HH34, where jets are observed 
with knots. Buhrke et al — explained that Kelvin-Helmholtz instability is the reason for knots in the 
jets. The jet must be ~ 10 times denser than its surrounding medium having velocity ~300 km/sec 
and Mach no. 30. Steady isolated jets may form structure through the growth of K-H modes. The 
stability properties of super magnetosonic astrophysical jets are subject of current interest. 

The linear theory of the combined effects of RT,KH and RM instabilities have been investigated 
earlier ^. Weakly nonlinear theoretical results of Kelvin-Helmholtz and Rayleigh- Taylor instability 
growth rates together with different aspects of density and shear velocity gradients have also been 
discussed E2- EH. 

In case of the temporal evolution of these instabilities nonlinear structures develop at the two fluid 
interface. The structure is called a bubble if the lighter fluid pushes across the unperturbed surface 
into the heavier fluid and a spike if the opposite takes place. The dynamics of such RTI and RMI 
generated nonlinear structures have been studied 031 — 033 under different physical situation using 
an expression near the tip of the bubble or spike up to second order in the transverse coordinate to 
unperturbed surface following Layzer's ^-approach. 

In the present paper, we investigate the combined effect of Rayleigh- Taylor, Richtmyer-Meshkov 
and Kelvin Helmholtz instabilities by extending the above method so as to include the effect of 
velocity shear induced contribution to the growth rate of the tip of the nonlinear mushroom like 
structures generated by shock wave (normal or oblique) incident on the unperturbed interface. 

In the event of excitation of RM instability due to normal incidence of shock in absence of velocity 



shear of the growth rate of the height of the finger like structures decay as It is however 

interesting to note that if the shock incidence is oblique (or if it passes across an irregular surface) the 
growth rate of the tip of the spiky structure does not decrease as (1/t) but attains a saturation value 
proportional to J k 2 (Av) 2 / (1 — At) where Av =difference is the tangential velocity of the fluids at the 
interface and At is the Atwood number. Thus the growth rate may be quite large if At — > 1 which 
may be likely in astrophysical situation and thus play an important role in formation of jets^,^. 

The paper is organized in the following manner. In section II is developed the basic equations 
describing the dynamics of nonlinear structures which evolve in consequence of the combined effects 
of these different types of hydrodynamical instabilities. In section III it is shown that the classi- 
cal results^^I] follow on linearization of the evolution equation describing the bubbles and spikes. 
Numerical as well as some analytical results regarding the saturation growth rates are presented in 
section IV. Finally section V presents a brief summary of this results. 

II. BASIC EQUATIONS OF EVOLUTION OF THE HY- 
DRO DYNAMIC INSTABILITIES 

Let the y — plane denote the unperturbed surface of separation of two fluids (the line y — in 
the two dimensional form of this problem). The fluid with density p a is assumed to overlie the fluid 
with density pb- The gravity g is assumed to act along the negative y- axis. Any perturbation of 
the horizontal interface or a shock driven impulse gives rise to Rayleigh- Taylor instability(p a > pb) 
or Richtmyer -Meshkov instability which in course of temporal evolution gives rise to nonlinear 
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interfacial structures. 

The two fluids separated by the horizontal boundary are further assumed to be in relative hori- 
zontal motion and thus subjected to Kelvin-Helmholtz instability arising due to horizontal velocity 
shear. Thus we are faced with the problem of the combined action of Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities. We shall see the same formulation will be applicable to Richtmyer-Meshkov 
instability associated with an oblique shock incident on the two fluid interface. 
After perturbation the finger shaped interface is assumed to take up a parabolic shape given by 

y(x, t) = r)(x, t) = r) (t) + r) 2 (t)(x - Vl (t)) 2 (1) 

For a bubble (here the lower fluid is pushing across the interface into the upper fluid with density 
p a > density pb) we have, 

r] > and r] 2 < (2) 

and for spike: 

r] < and r] 2 > (3) 

The height of the vertex of the parabola i.e, the height of the peak of the bubble (or spike) above 
the x-axis is |r/o(t)|. The position of the peak at time t is at x — i]i(t) and because of the relative 
streaming motion of the two fluids the peak moves parallel to the x-axis with velocity f]i{t). The 
densities of both fluids are uniform and fluid motion is supposed to be single mode potential flow. 

For the upper fluid with density p a we take the velocity potential 

4> a {x,y,t) = [a a (t) cos (k(x- Vl (t))) + f3 a (t) sin (k(x - Vl (t)))] e - k ^-^ - xu a (t); y>0 (4) 
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and for the lower fluid (density p b ) the velocity potential 

<f> h {x, y, t) = \a h (t) cos (k(x - Vl (t))) + b (t) sin (k(x - Vl (t)))} e fc (^°M) - xu b (t) + yb (t); y < (5) 

Before proceeding with the analysis of the kinematic and boundary conditions using the two fluid 
interface perturbation y = r](x, t) we forward the following justification for restricting the expansion 
to terms 0(x — r]i(t)) 2 . We are concerned only motion very close to the tip of the bubble or spike 
i.e., only in the region k \ x — rjx(t) |<< 1. Consequently one is justified in neglecting terms O 
(| x — f]i{t) |) 4 unless the coefficients of such terms are sufficiently large. Further it has been shown 
^231 that even it terms ~ T]±(t)(x — ^i(t)) 4 + r]e(t)(x — ^i(t)) 6 are retained the contribution from 
coefficients | 774 |, | T]q |<< that from | 7/ 2 | at least in the asymptotic state r — > 00. 

The kinematical boundary conditions satisfied at the interfacial surface y = 7](x,t)sae 
di] d(p a dr) d(j) a 



dt dx dx dy 

d(f) a dr] d(j) a d(fr b dr] 



(6) 



(7) 



dx dx dy dx dx dy 
The dynamical boundary conditions are next obtained from Bernoulli's equation for the two fluids 

+\{^<Paf +9ypa = -Pa + fa(t) (8) 



dt 2 

~ d ~tt + h^^ 2 + 9ypb = ~ Pb + m (9) 

by using the surface pressure equality 

Pa = Pb (10) 



leading to 



+ liy^a) 2 - k^fo) 2 + g( Pa - Pb )y 



dt dt ' 2 V Ta ' 2 

= fait) ~f h {t) (11) 



satisfied at the interface y = i](x,t) Now from Eq.(l) 

fto(t) - 2r h (t)r l2 (t)(x - m(t)) + 7) 2 (t)(x - ^{t)) 2 (12) 



dt 

Also utilizing the property that close to the tip of the bubble or spike, k\x — r)x(t)\ « 1, we express 
the velocity components in the following form 

Cl 1 

- = (««- kf3 a ) + k 2 a a (x - m) + f3 a k 2 ( V2 + k/2)(x - Vl ) 2 (13) 



dx 



v, 



dy 



- = ka a + k 2 (3 a (x - m) - k 2 a a (r] 2 + k/2)(x - r^) 2 (14) 



and similar expressions for v\, x and v b y 

Following Layzer's^ model we substitute for r] t , r] x , (v a ^) x , (v a (p)) y m the kinematic and boundary 
conditions represented by Eqs.(6),(7)and (ll)and equate coefficients of (x — r)i(t)) % ] (i = 0, 1,2) and 
neglect terms 0((x — t]i(t)) 1 ); (i > 3). This yields the following three algebraic equations for the three 
unknown b ,a b ,(3b : 

60 = " (3m-k/2) kaa (15) 

(377g + k/2) 

ab = w^m a (16) 

(r]2 + k/2)k/3 a -r] 2 (u a -u b ) 
Ph = H V2 - k/2) (17) 



and regarding the five other unknowns, viz r)o(t) , rj^t) , r) 2 (t) , a a (t) , /?&(£) the following five nonlinear 
ODE's [Eqs.(18)-(22)]. 



dr"" 6 



(18) 



(19) 



dr 



e 3 



A(6,r) (66-1) 



| 2(l-r)6(66~l) , N 2 fa,r) (66-1)6* 



+2 



L>i(6,0 A(6,r) 2^(26 -l) 2 

(46-i)(66-i) r, 



A(6,0(26-i) 5 



(K-H) 2 6-(K-H)(26 + 1)6 



(20) 



(26 - 1) 
£2(6, r) 



(/& - r /«) - r 



66 

26 



, 2(/ a - f b ) 

+ ^ — -?2 



^2(6,0 



(66 + 1)6 



2L> 2 (6,r)(66-l)(26-l) 



6 



4(K - H)(46 - 1) - f (28£ 2 2 - 46 - 1) 

6 



(21) 



dr a 



6(26 + 1) 
26 



(22) 



where 



6 = ki] ; 6 = V2/k; 6 = krjx 



(23) 



6 = k 2 a a / yfkg\ 6 = k 2 /3 a / yfkg, r = ty/~(kg) 



(24) 



V a = u a s[{k/g)] V b = u b yf{k/g)\ f a = / 



cir 



(25) 



The functions A?i i2 (6, r )) -^1,2(6) r ) where r = ^ is the density ratio are given by 



iVi(6, r) = 36(1 - r)i 2 2 + 12(4 + r)6 + (7 - r) 
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(26) 



A (6, r) = 12(r - 1)£ 2 2 + 4(r - 1)6 - (r + 1) 



(27) 



iV 2 (£ 2 , r) = 16(1 - r)£ 2 3 + 12(1 + r)£ 2 2 - (1 + r) 



(28) 



£> 2 (£ 2 ,r) = 2(l-r)£ 2 + (l + r) 



(29) 



The above set of five Eqs. (18)-(22) together with Eqs. (23)-(29)which define the different 
variables and functions describe the combined effect of RT and KH instabilities. 

On the other hand the impingement of an oblique shock on the two fluid interface causes the 
joint effect of Richtmyer-Meshkov and Kelvin-Helmholtz instability. The impact gives rise to an 
instantaneous acceleration which will change the normal velocity (y-component) by an amount Av = 
Vafter - Vbefore and transverse velocity (x-component) by Aw a(6 ) = (u a ^) after - (u a ( h )) before . Taking 
nonzero values only for the post shock velocities we replace the acceleration by their impulsive values. 
We set: 



and replace g — > AvS(t) The dynamical variables are non dimensionalized using normalization in 
terms of (kAv) instead of \fkg~. 

The combined effect of RM-KH instability resulting from oblique incidence of shock on the two 
fluid interface is then described by the same set of equations as Eqs.(18)-(22) together with the 
following replacements: 

(i) The second term on the RHS of Eq.(20)drops out. 

(ii) £3, £4 and r to be replaced by £3 = a a k 2 /(kAv), £ 4 = (3 a k 2 / (kAv) and r = t(kAv) respectively. 



dt 



u a S(t) -)■ Av(t) 



(30) 



10 



(iii) V a and V b by V a = u a /Av, V b = u b /Av. 

(iv) f a by 

7. = £=£-A W and A by A = ^ = ^A(T, (31) 
get At> or Aw 



III. LINEAR APPROXIMATION 

We now show that the usual combined RT and KH instability growth rates ^ are recovered on 
linearization of Eqs. (18)-(22). Let us put 



d(kr)i) d£ 5 ( p a){b) 

~ -r- = ot a V a + a b V b ; a a ^ {b) - 



or dr \ Pa + PbJ 

in Eq.(22)giving 

£ 4 = 2a 6 (K - H)ttAt ~ 2a ^ V - ~ V b)& ( 32 ) 

on linearization . In absence of velocity shear V a — V b = 0,we get £4 = 0. Thus the problem reduces to 

that of RT instability alone with no contribution from KH instability. Linearizing Eqs. (19),(20)and 
(21) we get 

f = (33) 



dr 



A T + a a a b {V a -V b f & (34) 



^ = -Pa(V a ~ V 6 )£ 3 (35) 
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At = P p~+p h is the Atwood number. Eq.(32) connecting £2 and £4 provides the consistency condition. 
The exponential growth rate due to combined effect of RT and KH instability coincides with the 
classical linear theory result ^2] 



(A) Combined effect of RT and KH instability: 

The growth rate of the RT instability induced nonlinear interfacial structures is further enhanced 
due to KH instability. Setting ^ = and ^ = the growth rate of the peak height of the bubbles 
and spikes are obtained by numerical integration of Eqs. (18)-(22) and the results are shown in Fig.l. 
The dependence of the growth rate on V a and V& keeping (V a — V ) unchanged are also indicated in the 
same diagrams.lt is found that for V b > V a the growth rate is greater than that for V a > V&(| V a — Vb | 
is the same for both cases); the asymptotic values is the two cases are however identical. Moreover 
for ^ > 1 Eqs.(18)-(22) show that as r — > 00 there occurs growth rate saturation given by 




(36) 



IV. RESULTS AND DISCUSSIONS 



ip \asym 



2A T 5(1 - A T ) 



(v a - v b y 



(37) 



\&)bubble 



\3(1 + A T ) 16(1 + Ax) 



and 



(6) 



asym 
spike 



2 At 5(1 + A T ) 



(V a ~ V b ) 2 



(38) 



\ 3(1 — At) 16(1- A t ) 



while 



(£0 



asym 
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for both bubble and spike respectively. 

Thus both saturation growth rate are enhanced for due to further destabilization caused by the 
velocity shear. 

On the other hand if r = p a /pb < 1(^4t = P p a +p b b < 0) there is no RT instability but it follows 
from Eqs.(37)and (38) that instability due to velocity shear (Kelvin -Helmholtz instability) persists 
on both the wind ward side and leeward side (i.e; both for bubbles and spikes) if (see Fig.2) 

and stabilized on both sides if (see Fig. 3 which shows oscillation of £iand £3 with respect to r) 

(V a -V b f<^(l-r) (40) 

If however (V a — Vb) 2 lies in the interval specified by the above inequalities, i.e; 

^(1 - r) < (V a - V h f < (r < 1) (41) 

it follows from the same two Eqs.(37)-(38) that the peak of the spike continues to steeper (instability) 
with r as the heavier fluid (density pb)pushes across the interface into the lighter fluid (density p a ) 
while the bubble height will execute low finite amplitude undulations. The above observation is 
shown to be suppressed in Fig. 4. At time t, the peak height which of the spike or the bubble occurs 
at x — ?7i(t) and thus moves to the right (x- increases) as i]i(r) increases with r. The spike peak 
height increase monotonically with t while that of the bubble undulates with low amplitude. The 
three dimension representation of the steepening of the peak of the spike as it moves along x-direction 
with time is shown in Fig. 5. In this respect there exists approximate qualitative agreement exists with 
the results of the weakly nonlinear analysis^^. 
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(B) Combined effect of Richtmyer-Meshkov and Kelvin-Helmholtz instability: oblique 
shock 

The time evolution of the two fluid interfacial structure resulting from the combined effect of 
Richtmyer-Meshkov and Kelvin-Helmholtz instabilities consequent to impingement of an shock is 
described by the set of Eqs.(18)-(22),(26)-(29) with modifications as shown in the set of Eq.(31). 
If the shock incidence is oblique then the normal component generates velocity shear and causes 
KH instability.^ The shock generated initial values of £ 3 and £ 4 are obtained from the impulsive 
accelerations represented by the 5— function terms in Eq.(30) giving 

"2(l-r)6(66-l ~ 



3)r=0 



(^4)r=0 



(42) 

(&)r=0 



(2£ 2 - l)(u b - ru a ) + 2£ 2 (w a - u b ) 



(43) 

J (&)r=0 



Av 

Results obtained from numerical solution of Eqs.(18)-(22) with modifications given by Eq.(31) 
subject to initial conditions (42) and (43) are presented in Fig. 6. The growth rate contributed 
in absence of velocity shear, i.e; by normally incident shock induced Richtmyer-Meshkov instability 
varies as t — > oo. However in presence of velocity shear the growth rate due to combined influence 
of RM and KH instability the growth rate approaches finite saturation value asymptotically. 

For RM-KH instability induced spikes it is given by the following closed expression 



fS .spike _ (dClypike 



5(1 + A T ) 
\ 16(1 -A T ) 



{u a -u b y/{Avy 



(44) 



which becomes large as the Atwood number A T — >•! (equivalentlyp a /p5 >> 1) 
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The following discussions suggest a higher plausibility of the effectiveness of the joint influence 
of RM and KH instability in the explanation of certain astrophysical phenomena. 

Corresponding to parameter values for Eagle Nebula(p a /p{, = 0.5 x 10 2 and \u a = 2 x 10 6 
cm sec -1 )^,^ the velocity of rise of the spike peak height h spifce (the height of the pillar) according 
to Eq.(44)is (4|)Soo ~ 0-^9 X 10 7 cm sec" 1 . Modification through inclusion of Rayleigh- Taylor 
instability effect (see Eq.(38)) can only slightly increases this value to ~ 10 7 cm sec -1 . This gives the 
time to reach the observed pillar height of 3 x 10 19 cm of the Eagle Nebula ~ 10 4 years. There are 
different time scales involved in the problem of development of the pillar of the Eagle Nebula. As 
pointed out by Pound, ^ there is a characteristic time scale for hydrodynamic motion Td yn ~ (A-u) -1 
where Av is the velocity shear inside the cloud. Corresponding to data given in ref.(3) this turns 
out to be Thydrodanmic ~ 10 5 yrs which is the upper time limit for development of the Eagle Nebula 
pillar ("elephant trunk"). But it is at least two orders of magnitude greater than the time scale 
Tcooi ~ 10 2 — 10 3 yrs imposed due to radiative cooling of the cloud E^ffl. j n comparison the time 
scale of the development of the pillar is found here ~ 10 4 yrs. Thus consequent to the hydrodynamic 
model based on the combined influence of Richtmyer-Meshkov and Kelvin-Helmholtz instability the 
gap between the two time scales t coo i and Thydrodynamicis reduced by one order of magnitude. 

A high Mach number, radiatevily cooled jet of astrophysical interest has been produced in lab- 
oratory using intense laser irradiation of a gold coneES. The evolution of the jet was imaged in 
emission and radiography. 

K-H instability growth rate has recently been observed in HED plasma experiment using Omega 
laser (A)=0. 351pm delivering 4.3 ±0.1kJ to the target overlapping 10 drive beams on to the ablator 
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^.Incompressible K-H growth rate peak to valley at Foam-Plastic interface has been compared with 
several analytical modes. 



V. Summary 

Finally we summarize the results: 

(a) If the heavier fluid overlies the lighter fluid the growth rate of both the bubble and spike peak 
heights due to RT instability are enhanced due to concurrent presence of velocity shear, i.e, K.H 
instability Fig.l. The asymptotic growth rates are given by Eqs. (37) and (38). 

(b) In the opposite case, i.e, if the overlying fluid is lighter and lower one is heavier (r = p a j p < 1) 
both the spike and bubble peak displacement increases continuously with time if 16 ^~ r - > < (V a — 
Vf>) 2 ,i.e, instability persists (Fig. 3) while stabilization occurs if (V a — V&) 2 < ( 16 (*~ r ) ) (Fig. 4 shows 
oscillation of peak heights of bubbles and spikes). 

(c) For 16( - ] L ~ r '^ < (V a — Vb) 2 < 16 ^~ r ^ for r< 1 the spike steepens with time (the peak height con- 
tinuously increases with time as indicated in Fig. 5. gives a three dimensional graph of displacement 
y against x and r). But the peak displacement of the bubble undulates within a small range Fig.4. 

(d) If the two fluid interface is subjected to an oblique shock Kelvin-Helmholtz instability due to 
generation of velocity shear occurs simultaneously with Richtmyer-Meshkov instability.The growth 
rates of bubbles and spikes due to this joint action are shown in Fig. 6. respectively. It is important 
to note that the growth rate of the combined action tends asymptotically to a saturation value given 
by Eq.(44); this is in contrast to that due to generation of RM instability due to normal shock 
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incidence for which the growth rate behaves as | as t— > oo. Moreover this growth rate as shown 
by Eq.(44) the rate of growth of this spike height has sufficiently large magnitude if the Atwood 
number A T — >l(p a « p&)-This may have interesting implication in the hydrodynamic explanation 
of formation of sufficiently long spiky jets in astrophysical situation, e.g, in case of the Eagle Nebula. 
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Figure 1: Initial values r = — = 1.5, £1 = —£2 = £3 = £4 = £5 = 0.1 for bubble;— £1 = 

Pb 

£2 = —£3 = £4 = £5 = 0.1 for spike. Plot showing variation of £1,^2, growth rate £3,^4 

and transverse displacement £5 of bubble and spike with V a = V& = 0.0 for solid black 

line-spike and broken black line for bubble. V a = 0.1, VJ, = 0.5 for broken blue line- 

20 

bubble and solid blue for spike, V a = 0.5, V& = 0.1, broken red line for bubble and solid 





Figure 2: r=^ = 0.4;Lower fluid denser. Dashed line for spike (heavier fluid pushes 
into lighter fluid) and unbroken line for bubble. V a = 0.8, V b = —0.6. Initial condition 
as in Fig.l. 

and for following relation < (V a ~ Vb^ 




Figure 3: r=OA,V a = 0.0, H = 0.2.(14 - V b f < ^jf 1 .Initial condition as in Fig.l. 
Unbroken line for bubble and dashed line for spike. 



22 



0.1 



0.05 



-0.05 



-0.1 



-0.2 
20 



0.5 - 



* 0.2 



-0.2 





Figure 4: r=0A,V a = 0.6, H = -0.4; < (V a - V b f < .Initial condition 

as before (Fig.3.). Unbroken line for bubble and dashed line for spike;height of spike 
peak increases monotonically with time (steepening) ;bubble depth undulates. 
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Figure 5: 3 dimensional plot of spike (Interface Y = ?7o(t)+^2(t)(^ — ?7i(T)) 2 )belonging 
to the plot given in fig. 4. 
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Figure 6: Oblique shock: RM and KH instability for spike (dashed line) and bubble 
(unbroken line). Initial values as in fig.l. and V a — 0.1, V b — 0.5. 
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